Importance of the Wick rotation on Tunnelling 
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A continuous complex rotation of time t te~ l9 is shown to smooth out the huge fluctuations 
that characterise chaotic tunnelling. This is illustrated in the kicked rotor model (quantum standard 
map) where the period of the map is complexified: the associated chaotic classical dynamics, if 
significant for 9 = 0, is blurred out long before the Wick rotation is completed (9 = 7r/2). The 
influence of resonances on tunnelling rates weakens exponentially as 9 increases from zero, all the 
more rapidly the sharper the fluctuations. The long range fluctuations can therefore be identified 
in a deterministic way without ambiguity. When the last ones have been washed out, tunnelling 
recovers the (quasi-)integrable exponential behaviour governed by the action of a regular instanton. 
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It is often admitted that the Wick rotation t — > —it 
provides a straightforward route from a Minkowskian 
metric to an Euclidean metric or from a zero temper- 
ature quantum model to a finite temperature statistical 
model. In this note we will show that even a small com- 
plex rotation of time, t — > e~ w t can affect drastically the 
spectral properties of some operators that encapsulate 
the quantum dynamics. In particular, the fluctuations of 
tunnelling rates and therefore some quantum transport 
properties, can change by several order of magnitude even 
if 9 < 1. 

This work is especially motivated by the topic of 
chaotic tunnelling, a theoretically challenging subject 
that has been resisting a systematic and satisfying so- 
lution for more than sixteen years. Since the pioneer 
quantitative studies [l|, it has been widely observed @ 
that tunnelling, i.e. a quantum process that is forbid- 
den at a classical level, exhibits a huge sensitivity to 
any perturbation and generically fluctuates by several or- 
ders of magnitudes if the underlying classical dynamics 
is non-integrable. Understanding, predicting and con- 
trolling this behaviour, remain widely open problems for 
which several strategies have been proposed. Following 
semiclassical methods, known to be successful in tack- 
ling the issues of tunnelling in integrable systems, the 
most natural strategy is to try to express a tunnelling 
rate (for instance a decay time or an oscillation period 
between two symmetric wells) in terms of complex so- 
lutions of Hamilton's equations only. Already highly 
non-trivial in the multidimensional (non-separable) in- 
tegrable or quasi- integrable cases 0, @, this program 
appears even more difficult when some of the integrals 
of motions are strongly broken. Indeed, it has been 
discovered 0] that chaos reveals itself in the complex 
phase-space though some fractal structures, the so-called 
Laputa islands, that look like agglomerates of complex 



classical trajectories. It is only recently that some en- 
couraging, significant steps were realised in retaining the 
relevant semiclassical skeleton Q for tunnelling. To by- 
pass this purely semiclassical strategy, it has been pro- 
posed [H to replace the chaotic, though deterministic, 
transition amplitudes (i.e. involving one or more chaotic 
quantum states) by some random matrix elements and 
hopefully compute an average tunnelling behaviour on 
the appropriate statistical ensemble. This hybrid ap- 
proach, mixing together some semiclassical integrable in- 
gredients and statistical chaotic ones, allows to establish 
a precise and quantitative connection between the quan- 
tum and classical resonances which play a crucial role for 
understanding multidimensional tunnelling When 
classical chaos is well developed, the overlap of classi- 
cal resonances has a quantum counterpart: the collective 
effect of the quasi-coincidences of quantum frequencies, 
which cannot be isolated one from the other, not only 
causes the huge fluctuations of a tunnel ling rate but also 
enhances its average behaviour [Til, [U, EDf • In this 

con- 
text, the complex rotation of time that is proposed in this 
note provides a continuous way to understand how the 
resonances conspire to create a chaotic tunnelling regime. 
No statistical ingredients are required and therefore the 
present approach establishes a new bridge between the 
two strategies described above since it provides a deter- 
ministic, resonances-governed transition between a reg- 
ular and a chaotic tunnelling regime where the complex 
classical solutions play a natural and essential role. 

One of the simplest models where chaotic tunnelling 
is at work corresponds to a quantum system with one 
degree of freedom whose dynamics is a sequence of pe- 
riodically alternating kinetic and potential motions [14j . 
After one period r, the quantum evolution operator (also 
known as the quantum map or the Floquet operator) is 
given by 



U(r) 



e -iTf(p)/h e -iTg(q)/h 



(1) 
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The associated classical dynamics after one period is de- 
scribed by the discrete Hamilton equations (the Poincare 
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map): (po,go) ^ wit h Pi = Po - T~g'(q ) and 

<Zi = <Zo + T f'(Pi) (the primes stand for the derivative of 
the smooth functions / and g). The choice f(p)=p 2 /2 
and (7(g) =7 cos q where 7 is a real parameter, corre- 
sponds to the well known kicked rotor model (the stan- 
dard map) that is extensively studied for understanding 
the subtle interplay between classical and quantum trans- 
port [151 ]. It is convenient to work with the usual sym- 
metric double-well situation that can be obtained from 
the kicked rotor model by unfolding the cylindrical phase- 
space on a double spatial period and work in the following 
with g(q) = 7 cos 2q with strictly periodic quantum peri- 
odicity for any g-translation by 2n [l(| . The classical dy- 
namics is shown in Figurc[T] No real trajectories connect 




FIG. 1: Real phase space trajectories for the kicked rotor 
with r real and 7 = 0.25. The upper panel corresponds to 
the quasi-integrable case (r — > + ). The (red) thick line is 
the trajectory joining the unstable fixed point (p = 0, q = 0) 
to (p = 0, q = it) after an infinite number of iterations of the 
map. The mixed dynamics in the lower panel corresponds 
to r = 1 where a chaotic sea separates the two regular islands 
in the neighbourhood of the stable points (p — 0, q — ±7r/2). 

the interior of the two symmetric islands in the neigh- 
bourhood of the two fixed points p — and q — ±7r/2 
that remain stable as long as 

r e ] -1/V7, W7[. (2) 

At the quantum level, we will consider the dynamics that 
is described in terms of the eigenstates of (TT|) that are 
strictly invariant under a g-translation by 27r. The classi- 
cal phase-space parity symmetry (p, q) 1— > (— p, —q) allows 
to classify the spectrum of ([1]) accordingly: we will de- 
note it+ (resp. u~) the eigenvalues corresponding to the 
symmetric (resp. antisymmetric) eigenstates \4>n)i n is 
an quantum number labelling the doublet. For real r, the 
quasi-energies defined modulo 2-kK/t by = ^ hi(w^) 
are real. For a given n, when \<j)^) have their Husimi 



representation localised inside the two wells, tunnelling 
is characterised by the splitting Ae n = e~ — e+. In the 
following, we will drop the n, keeping in mind that we 
will work with the central doublet i.e., selected by the 
criterion of having the maximal overlap with a coher- 
ent state that is located on the stable fixed points p7| . 
The expected semiclassical behaviour of Ae in a (quasi- 
)integrable regime where no classical resonance can be 
resolved by quantum eyes is given by Q 

Ae ~ ah»e- lAl/h , (3) 

where a and A are made of classical (ft-independent) in- 
gredients. The exponent v depends on the nature (in- 
tegrable or quasi-integrable) of the dynamics. For a Id 
time-independent system it is well known [l8j that v = 1 
and both a and \A\ can be interpreted in terms of instan- 
tons [19(, i.e. expressed from the real trajectory joining 
(in an infinite time) the two fixed points once the Wick 
rotation (9 = tt/2) has been completed. In multidimen- 
sional systems or in the case of a quantum map, the 
generic presence of quantum (and classical) resonances 
will drastically modify ([3]). This can be understood as 
follows: for r real, the it's accumulate on the unit cir- 
cle and two quasi-energies differing by almost an integer 
multiple of 2ttH/t will make the corresponding u's almost 
coincide and generate small denominators in any pertur- 
bative expansion. Some transitions between the central 
doublet and some excited doublets are therefore greatly 
enhanced and give birth to large fluctuations when a con- 
trol parameter is varied. In figure [31 the latter is chosen 
to be the inverse of the (effective) Planck constant and Ae 
is plotted in black dots for r = 1 and 7 = 0.25 corre- 
sponding to the classical dynamics in the lower panel in 
figure [TJ 

The operator U(r) can be analytically continued in 
the lower half-plane of complex r's. Giving a negative 
imaginary part to r makes U non-unitary and its eigen- 
values escape from the unit circle into the whole complex 
plane. In the quasi-integrable case obtained for small 
enough |r| and 7 for the two exponentials in |T]) to almost 
commute, it is expected that the eigenvalues it„ remain 
near the logarithmic spiral e~ 1TS / h parametrised by the 
real s. It is surprisingly the case even in a regime where 
classical chaos is well developed (for |t| = 1, 7 = 0.25, 
see figure [IJ and r being far from the real axis. A 
complex rotation of r by a positive angle 9= — argr 
appears to be a simple tool to unfold the unit circle 
into a spiral and therefore to move the resonant dou- 
blets away from each other: an increase of the quasi- 
energy by 2ttH/\t\ is accompanied by a shrinking of the 
modulus of the it's by a factor of exp (— 27rsin6'). The 
larger 9, the less the influence of the resonances on the 
tunnelling splitting. This is confirmed by the plots in fig- 
ure [3J where we can see how the fluctuations of |Ae| are 
progressively smoothed by a complex time rotation. For 
6*'s about 10~ 3 , only the more acute spikes are eroded. 
For 9 ~ 0.01-7T, only the long range fluctuations in l/h 
have survived and the staircase-like structure that can 
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FIG. 2: (Color online) Eigenvalues w„ 
plane for 7 = 0.25, ft ~ 1/8.002 and, 



of fTJ in the complex 
t = e~ ie with 9 = 



large dots (black) ; 9 = 10 71": small dots (green); 9 = tt/10 
stars (blue); 9 = n/4: diamonds (magenta); 9 = 0.457T 
squares (cyan). The tunnelling central doublet is circled (red). 
Its splitting is of order at most 10 -5 (see figure [5} and there- 
fore cannot be resolved at this scale. The continuous line are 
the logarithmic spirals e~ 1TS ^ h parametrised by the real s. 




15 20 25 
2%/h 



FIG. 3: (Color online) Modulus of the central tunnelling split- 
ting for the kicked rotor with 7 = 0.25, r = e _lfl . 9 = 0: dots 
(black) ; 9 = 10" 3 : thick solid line (cyan) ; 9 = 2 10~ 3 : thin 
solid line (magenta); 9 = 3 10 -3 : starred solid line (blue); 
9 = 3 10 -3 7r: dotted line (green); 9 = 10 _2 7r: dashed line 
(yellow); 9 — 5 10~ 2 7r: squared line (violet); 9 = 7r/10: dia- 
mond line (orange); n/5 ^ 9 ^ n/2: thick solid line (red). 



be observed on the semi-logarithmic plot is very similar 
to the average behaviour of In | Ae| that is obtained from 
the resonant-assisted tunnelling hybrid method described 
above (compare with [ll], Figures 1 and 2] and [H, Fig- 
ure 10]) but here no averaging process is required and 
there is no rough transition between staircase steps due 
to the semiclassical truncation of the Hamiltonian be- 
tween a regular part and the random matrix that mod- 
els the transitions involving a chaotic state. What is 
unexpected and of course requires further future inves- 
tigations, is that the splitting curves (in semi-log plots) 
converge quickly, say for tt/5 < 9 < tt/2, to a straight lino 
in agreement with the (quasi-)integrable behaviour ([3]). 
The slope — 1^4| is roughly the same as the one that ap- 
pears in between the plateaus for smaller 9. Moreover, 
keeping |r| = 1 and computing the slope of the limiting 
straight lines obtained for several values of 7, including 
the ones for which the stable islands are almost dissolved 
in the chaotic sea for real r (according to Q the bifurca- 
tion point where the stable fixed points lose their stability 
occurs at 7 = 1). It is shown in figure[4]that |^4| ~ 4^/% 
which can be interpreted as the action of the instanton 
of the regular regime. It is easy to check that the com- 
plete Wick rotation (9 = tt/2) of the kicked rotor corre- 
sponds to a simple g-translation by tt/2 of the dynamics 
obtained for 6 = 0. Therefore the action of the instanton 
that joins (p = 0, q = ±tt) in an infinite purely imaginary 
time is the area 4.^/7 under the separatrix shown in the 
upper panel of figurc[T] From a purely semiclassical point 
of view, as recalled above, one needs to elucidate the rel- 
evant structures in the complex phase-space. For 8 > 0, 
those are expected to be more tractable, compared to the 




FIG. 4: (Color online) As 9 increases, the semi-logarithmic 
graphs of 1/h 1— > [Ae| eventually accumulate on straight lines 
(shown in the inset) whose slope —\A\ (sec Q) is reported 
for several values of 7 (|r| = 1 being fixed). The red arrows 
indicate the case 7 = 0.25 corresponding to figures [1] and [3] 
The continuous black line is the graph of —4^/7 

situation where 9 = 0. It can be predicted without too 
much risk that the inextricable fractal structure of the 
Laputa islands (for the kicked rotor see Figure 2 of Q) 
simplifies when increasing 9. In other words, decreasing 9 
towards should allow us to understand the formation 
of the delicate pattern of the Laputa chains from isolated 
orbits in the same time as we see in figure |3] how the re- 
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sonances overlap and glue together in a dense set of Ae 
fluctuations. 

The above analysis raises of course several (still open) 
questions. In particular the role of the regular instanton 
is puzzling since it represents a primitive classical trajec- 
tory with infinite period and therefore seems out of reach, 
if one were to tentatively try to obtain a trace formula a 
la Gutzwiller starting from the estimate 



Ae 



2ih tr(SU N {T)) 
tr(U N (r)) 



Nt 



J2 (^.-iNTtt/h _|_ e -iNTe~ /f^J 



(4) 



(5 is the parity operator) valid in a regime where the 
integer N is sufficiently large for the central doublet 
to dominate all the others but sufficiently small so 
that NrAe/h <C 1. Indeed, the usual resummation tech- 
niques lead to an expansion arranged in ascendingorder 
of the length of the primitive (pseudo-)periodic [211 ]. As 



far as tunnelling is concerned, I have not been able to ex- 
tract any relevant information from the shortest complex 
pseudo-periodic orbits. The chaotic character that seems 
to fade away for < 9 < ir/2, must also be understood 
better. In this interval, all the orbits that are stable 
for 8 = or ir/2 get generically destabilised since the 
traces of the monodromy matrices escape from ] — 2,2[ 
into the complex plane but the mixing properties in a 
non-compact complex phase space may vanish as well. 
The non vanishing Imr has the flavour of a dissipation 
and the suppression of the chaos-enhancement of tun- 
nelling reminds of the Caldeira-Leggett model [2(| ■ 

The continuous complex rotation of time undoubtedly 
opens new promising perspectives for gaining a proper 
understanding of chaotic tunnelling and, as recalled in 
the introduction, this new approach sheds light on the 
subtle non-analytical connection that may happens con- 
cerning the transport properties of complex systems that 
are linked by the Wick rotation. 

I thank O. Brodier, D. Delande, B. Gremaud for con- 
tinuous stimulating discussions on this subject and S. 
Nicolis for his careful reading of the manuscript. 
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